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We present a detailed study of the population imbalanced three-component Hubbard chain with attractive 
interactions. Such a system can be realized experimentally with three different hyperfine states of ultra cold 
^Li atoms in an optical lattice. We find that there are different phases that compete with each other in this 
system: A molecular superfluid phase in which the three fermion species pair up to form molecules (trions), 
a usual pairing phase involving two species with exactly opposite momenta, and a more exotic generalized 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase consisting of three competing pairing tendencies with differ- 
ent non-zero center-of-mass momenta. At large attractive interactions the system exhibits strong tendencies 
towards collapse and phase separation. Employing the density-matrix-renormalization-group-method (DMRG) 
to determine the decay exponents of the various correlators we establish the phase diagram of this model for 
different fillings and interactions. We also discuss the experimentally relevant situation in a trap and report the 
existence of an additional region where two species are dynamically balanced. 

PACS numbers: 03.75.Ss, 03.75.Mn, 71.10.Fd, 71.10.Pm 



Introduction - Pair formation of fermions is an impoitant 
phenomenon occurring in a variety of physical systems rang- 
ing from superconductors and superfluids to dense quark mat- 
ter. The pairing of two balanced species of fermions is well 
understood within the Bardeen-Cooper-Schiieffer (BCS) the- 
ory [1], but pairing and molecule formation involving three 
flavors are still intriguing, especially if the Fermi surfaces 
of the different flavors do not match. Recent advances in 
methods for trapping and controlling ultracold atoms have 
opened up the possibility of experimentally observing a three- 
component mixture of ultracold ^Li fermionic atoms with 
attractive interactions [2, 3]. These experimental prospects 
have fueled a variety of theoretical studies of systems with 
equal densities of flavors, in which the formation of three- 
body molecules is favored for strong attraction. In one dimen- 
sion, the molecular superfluid prevails for all attractive inter- 
actions [4], whereas in higher dimensions, an additional color 
superfluid phase has been found for smaller attractions [5]. 

In the more general, but still experimentally realizable case, 
where the flavour densities are different, it is an open question 
how the system maintains its binding tendencies and whether 
it favors trion or pair formation. In the two-flavor case - as 
exemplified by the imbalanced attractive Hubbard model - 
an extended Fulde-Ferrell-Ovchinnikov-Larkin (FFLO) phase 
with a well defined momentum Q = {kp^ — ^fxI of the 
Cooper pairs has been found [6]. In the three-flavor case there 
are three different possible FFLO wave vectors and it is un- 
clear which one will dominate the pairing properties. Moti- 
vated by these fundamental theoretical question and their rele- 
vance to upcoming experiments, we have investigated the case 
of generic filling in the three-component Hubbard chain with 
attractive interactions - the simplest possible model describing 
these binding tendencies. In this Letter, we present the pairing 
phase diagram for different imbalance and interactions. 

Model - We concern ourselves with a population imbal- 
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FIG. 1: (Color online) Density configurations for a given filling n 
form a triangle in ni-n2-n:i space. Because of symmetry, it is suffi- 
cient to consider only the subset ni < n2 < 713, shown on the right 
hand side. Dots represent the systems studied in this work and the 
three encircled numbers refer to selected configurations presented in 
Figs. 2 and 4. 



anced three-component Fermi gas in a one-dimensional (ID) 
optical lattice, which can be described by a Hubbard Hamil- 
tonian 
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1, 2, 3 on site / and ni,a - 
cupation number operator. The parameters t, which we set to 
t 1, and Ua/s = U < characterize the nearest-neighbor 
hopping and the SU{3) invariant attractive on-site interac- 
tion, respectively. The flavor density configurations compati- 
ble with a given total density n :— 1/L '^Q' where L is the 
length of the chain, form a triangle in ni-n2-n3-space. Due 
to permutation symmetry, it is sufficient to peiform calcula- 
tions on the subset rii < n2 < n^, see Fig. 1. The centroid is 
the point where all densities are equal. It has been shown [4] 
that any U < leads to the formation of a Luttinger liquid 
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TABLE I: Important correlation functions in the three-component 
attractive Hubbard chain. 



Fluctuation 



Expression 



Gail, I') Green's functions 
Paf3{l,l') a/3-pairing correlations 
T{l,l') trion correlations 



of three-flavor molecules (trions). Along the medians, two 
out of the three densities are equal. Depending on whether 
7ii < = 713 or Til = < "-3 holds. One encounters 
either a situation of flavor- 1 fermions mixed with flavor (2- 
3) cooper pairs (a Bose-Fermi mixture) or trions mixed with 
flavor-3 fermions (a Fermi-Fermi mixture). Along the edges 
of the triangle one recovers the physics of the two-flavor case. 

General considerations - In this ID system , long-ranged 
pairing correlations are suppressed by phase fluctuations, 
leading to an algebraic or exponential decay. In the case 
of algebraic decay, the leading asymptotic behavior of the 
equal-time correlators is described by power laws C(r) c>c 
cos (fcr) /r'^ , where k, an appropriate integer combination of 
the Fermi momenta fcj? „, is the ordering wave vector, and all 
the parameters depend on the type of fluctuation considered. 
The physics is determined by the dominant correlation with 
the slowest decay and thus the smallest exponent v. The rel- 
evant correlations for attractive interactions considered in this 
Letter are summarized in Tab. I. 

Correlations - To determine the dominant fluctuations, we 
use the density-matrix-renormalization-group (DMRG) algo- 
rithm [7, 8] to calculate the correlation functions indicated in 
Tab. I for open chains with L — 128 sites at total density 
71 = 3/8 and 3/4. The correlation exponents are subsequently 
extracted by fitting the power law decays of the average cor- 
relators C'(r) — C{l,l + r)/{L ~ r), to reduce boundary 
effects. 

Two representative examples for density 77 = 3/4 and in- 
teraction U/t ~ —4 are shown in Fig. 2. In the situation 
depicted in the left column, two species are balanced, while 
the remaining flavor remains a minority (771 < 712 = 773). 
This corresponds to a two-component (Bose-Fermi) Luttinger 
liquid with expected central charge c— 2, see Fig. 1. The bal- 
anced species form BCS pairs with zero momentum, as can 
be deduced from the strong peak at A: = in the Fourier 
transform of the pairing correlators shown in the inset. Ac- 
cordingly, single particle excitations in these two channels are 
gapped and the corresponding Green's functions decay expo- 
nentially. Note that the momentum of the pairing correlations 
of the (2-3) pairs is not affected by the presence of the specta- 
tor flavor (1). The correlation exponent is however increased 
in comparison to the two-component model. Although trionic 
correlations also decay algebraically (here with wave vector 
kp i), they decay significantly faster than the (2-3) pairing 




FIG. 2: (Color online) Absolute values of Green's functions and 
pairing correlations for two systems at density n = 3/4 and inter- 
action U/t = —4. The left column shows a balanced-imbalanced 
situation with TS! = (A'^i, A^2, A^s) = (16,40,40) particles, while 
the right column illustrates the generic case with N — (12, 36, 48) 
fermions. The encircles numbers refer to Fig. 1. The inset shows the 
Fourier transform of the real space correlators presented in the main 
panels. The exponents are extracted from power law fits. 



correlations, which thus govern the physics of this regime. 

If all three densities are different, the system behaves as 
a three-channel Luttinger liquid with expected central charge 
c = 3 in which all Green's functions and pairing correlators 
decay as power laws. An example of such a system is shown 
in the right column of Fig. 2. Similarly to the imbalanced two- 
component Hubbard model [6], the system now forms Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) pairs with finite center- 
of-mass momentum. The important difference however is 
that here, there are three possible FFLO pairs with momenta 
Qap _ _ kp^pl and correlation exponents v"'^. This 

is nicely illustrated by the distinct peaks at fc = Q"'' 7^ 
in the Fourier transform of the three pairing correlations. In 
this example, the pairing in channel (2-3) is dominant, pairs 
involving species (1-2) are subdominant, and the (1-3) pairing 
channel decays very quickly. This sequence suggests an in- 
crease of the correlation exponent with increasing pairing mo- 
mentum. Similar to the Bose-Fermi case discussed before, the 
molecular superfluid correlations decay more quickly than the 
dominant FFLO fluctuations, but they appear to decay more 
slowly than the second pairing channel. Multiple FFLO wave 
vectors have recently been reported also for two-chain sys- 
tems [9]. 

Phase diagram - The complete numerical pairing phase di- 
agram for total density 77 = 3/4 and weak to moderate interac- 
tions U/t < 4 is presented in Fig. 3. Each row (a)-(c) contains 
a small triangle highlighting the dominant correlation func- 
tion, while the larger triangles quantify the dominance based 
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FIG. 3: (Color online) Dominant correlations of the three-component Hubbard chain at filling n = 3/4 for various interactions U. The 
smaller triangles indicate the leading correlation function i.e. the one with the smallest decay exponent, while the larger triangles contain an 
additional color coding of the contrast and reliability of this determination, see text. The more intense the color, the smaller the overlap with 
the subdominant exponent and the more reliable the prediction. Row (a) is restricted to the decay of the three Green 's functions, in row (b), the 
dominant pairing is indicated, while row (c) represents the combined phase diagram, including both Green's function and pairing correlation 
exponents. 



on the fitted exponents and their uncertainty [16]. Let us first 
focus on the Green's functions shown in row (a) of Fig. 3. For 
small interactions, our results clearly reveal that the Green's 
functions of the majority species decay most slowly. With 
increasing U, one observes the growth of an additional re- 
gion centered around the lower part of the median (region A 
in Fig. 3, indicated in the upper rightmost triangle). In this 
area, most of the majority species are paired up due to strong 
attraction. Since the density of the remaining free fermions 
is smaller than the density of the nominal minority flavor, the 
Green's function of the latter one has the slowest power law 
decay. In row (b) of Fig. 3, we present the dominant pairing 
correlations. The formation of trions is only favored at the 
SU(3) symmetric centroid. For generic density configurations 
and weak interactions, we find dominant FFLO pairs formed 
by the minority species, except along the medians, where BCS 
{Q = 0) pairing is clearly dominant. In the weak coupling 
limit \U\ 0, we thus expect the qualitative phase diagram 
sketched in Fig. 3[row (b), column (i)]. In this limit, it is 
energetically favorable to have as many free fermions as pos- 
sible, hence only the minority species form pairs. For stronger 
interactions, the patches around the bases of the medians (re- 



gion B) spread out, hereby reducing the areas at the top (re- 
gion C). Based on this observation, one would conjecture that 
in the strong coupling limit \U\ oo, majority pairs domi- 
nate, as sketched schematically in Fig. 3 [(b), (iv)]. However, 
based on a perturbative argument mt/U, one realizes that the 
hopping of trions, pairs, and individual fermions intervene at 
third, second, and first order, respectively. This tunable hop- 
ping asymmetry leads to different sorts of collapses, demix- 
ing, and phase separation at large \U\/t, reminiscent of the 
hopping asymmetric two flavor Hubbard model [10]. In ad- 
dition, we find evidence for the formation of multi-particle 
composite liquids close to phase separation, as advocated re- 
cently in [11]. At intermediate couplings, the boundaries of 
the dominant correlations are surprisingly well described by 
lines of equal density difference, as sketched in Figs. 3 [(b), (ii) 
and (iii)]. Trionic fluctuations are never dominant away from 
the centroid, but their correlation exponent decreases substan- 
tially with increasing attraction. 

Dimensional crossover - Following the standard analysis 
of the dimensional crossover, see, e.g., Ref. 12, we compare 
the dominant two-particle correlations with the single parti- 
cle Green's functions and argue that upon lowering the tem- 
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FIG. 4: f Color online j Local densities for the three-component Hub- 
bard chain at filling n = 3/8, with {Ni, N2, N3) = (12, 16, 20) i« a 
Trap 'Krap ~ 0.002 for various interactions U. The encircled number 
refers to the density configuration indicated in Fig. 1. 

perature of the system, one either enters a phase dominated 
by the physics of the single chain FFLO correlations if the 
two-particle correlations decay more slowly than the Green's 
function, or, one first flows to a higher dimensional Fermi liq- 
uid if the opposite condition is satisfied. Combining rows (a) 
and (b) leads to the complete phase diagram shown in row (c). 
For weak attraction, the fermion Green's functions are dom- 
inant, except for the centroid and the midpoint of the edges. 
Pairing correlations become more important with increasing 
U, first only along the medians, where fc = Cooper pairs 
are formed, but eventually, one also observes FFLO pairing 
around the center of the triangle. From this sequence of im- 
ages, we conclude that exotic pairing tendencies carry over to 
the quasi- ID case in situations where all three species have 
similar densities. We have also undertaken the same analysis 
for total density n = 3/8 (not shown) and found a very simi- 
lar behavior, albeit with a crossover from the Fermi liquid to 
the pairing phase taking place at weaker interactions. 

Influence of the trap - In experiments, the atomic cloud is 
confined in a harmonic trap of the form 

H,..v^vY,(l-^^\ n,. (2) 
1=1 ^ ^ 

with a trapping potential V. For the case of three equal flavor 
densities an atomic density wave has been reported in [13]. 
For the flavor imbalanced case an interesting effect of the trap 
is illustrated in Fig. 4, for a particular filling and realistic val- 
ues of U and V . For small U, the area around the center of 
the trap resembles the corresponding homogeneous system. 
We thus expect to find the multiple FFLO behavior discussed 
for the homogeneous system also in the presence of a trap. 
However, for stronger interactions, the system rearranges the 
particles in such a way that two flavors have equal densities at 
the center of the trap, hereby forming polarized wings. The 
pairing behavior remains, but the FFLO signatures disappear 
and one finds BCS pairs with zero momentum instead. The 
phase diagram presented in Fig. 3 thus becomes even richer 
after the inclusion of the trap, exhibiting an additional BCS 
dominated strong coupling regime. At higher filling, the trap 
can also lead to the formation of inert inner shells, in which 
the density of the majority flavor is exactly equal to one. A 
similar effect has been observed in the two-component Hub- 



bard chain [6]. 

Conclusion - We have established the pairing phase dia- 
gram of the three-component Hubbard chain with attractive 
interactions using the density-matrix-renormalization-group 
algorithm. At generic imbalance there is a coexistence of mul- 
tiple FFLO modes. It will be interesting to investigate the fate 
of this coexistence in higher dimensions, where Ward iden- 
tities predict phase separated mixtures of paired states [14]. 
The phase diagram of the present non-integrable lattice model 
is richer than the one reported recently for an integrable three 
flavor continuum model [15], due to the formation of multi- 
particle composites, collapse and phase separation at strong 
attraction. We are presently investigating the precise phase 
boundaries of these density instabilities. 
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